Correlation Lengths in Quantum Spin Ladders 
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Analytic expressions for the correlation length temperature dependences are given for antiferro- 
magnetic spin-1/2 Heisenberg ladders using a finite-size non-linear cr-model approach. These cal- 
culations rely on identifying three successive crossover regimes as a function of temperature. In each 
of these regimes, precise and controlled approximations are formulated. The analytical results are 
found to be in excellent agreement with Monte Carlo simulations for the Heisenberg Hamiltonian. 
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Low-dimensional quantum magnets exhibit many 
novel collective low-temperature properties. An inter- 
esting recent development is the system of spin ladders 
that are arrays of n coupled antiferromagnetic Heisen- 
berg chains. Experiments and numerical studies have re- 
vealed many fascinating aspects of these systems, such as 
the surprise that spin-1/2 ladders composed of an even 
number of chains have a gap in the excitation spectrum, 
while those with an odd number of chains are gapless Q . 

For the two-dimensional square-lattice spin-1/2 
Heisenberg antiferromagnet with nearest-neighbor inter- 
actions, theoretical predictions for the spin-spin cor- 
relation length based on the (2+l)-dimensional non- 
linear cr-model |S| are in excellent quantitative agree- 
ment with neutron scattering results ||] . In Ref. || , one 
of us pointed out that it is useful to view a spin lad- 
der as a finite sized two-dimensional antiferromagnet, 
instead of a system of coupled spin chains. The zero- 
temperature properties could then be readily obtained 
from the finite-size scaling of the (2+l)-dimensional 
non-linear cr-model, where one of the spatial directions 
is finite and the remaining directions infinite, one of 
them being the Euclidean time direction. In this manner 
one could obtain both a clear picture of the dimensional 
crossover to the two-dimensional square-lattice antifer- 
romagnet, as well as precise analytical estimates of the 
zero-temperature correlation lengths and gaps in spin 
ladders. Moreover, it was predicted that the spin gap, 
A, is simply related to the spin-spin correlation length, 
£, such that 
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where J is the antiferromagnetic coupling, a is the lattice 
spacing, and c is the physical spin wave velocity of the 
two-dimensional square lattice antiferromagnet at zero 
temperature. Applied to S — 1/2, the right hand side of 
Eq. (0) should be 1.68. 

At first sight, an analysis at non-zero tempera- 
tures similar to that at zero temperature may appear 
formidable. One of the purposes of this paper is to show 
that this is not so. Once one recognizes three distinct 



crossover scales as a function of temperature, the temper- 
ature dependence of the correlation length can be com- 
puted with equal ease. In the language of the cr-model, 
the system is a box whose two sides L y = na (the spatial 
direction along the width of the ladders) and (3hc (the 
Euclidean time direction) are finite in extent, where j3 is 
the inverse temperature. Successive crossovers take place 
as the system goes through regimes where (3hc 3> L y , 
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L y , and f3fic 



(There is an additional 



crossover at even higher temperatures, where the system 
becomes entirely classical, which is not discussed here.) 
In each of these three regimes, it is possible to formu- 
late precise and controlled analytical methods to com- 
pute correlation lengths of spin ladders. These estimates 
are then tested against new numerical simulations, and 
we find excellent agreement between the numerical and 
the analytical results. Because we use periodic boundary 
conditions across the width of the ladders in the cr-model 
approach ||, it is necessary for us to carry out simula- 
tions with periodic boundary conditions as well. How- 
ever, we show that the choice of boundary conditions 
hardly matters for ladders of width larger than n = 4. 
Finally, we show that the zero-temperature predictions 
of the correlation lengths that contain no adjustable 
parameters, are in good agreement with the simulation 
results, and that the relation shown in Eq. (0) holds with 
remarkable accuracy. Thus, we believe that we have a 
reasonably complete picture of spin fluctuations in spin 
ladders. 

The Hamiltonian for a Heisenberg ladder is 



H = J\ 



<«j> ii 



(2) 



where S; = ^cr^ is the quantum spin operator at site i, 
while and (ij)± denote nearest neighbors along and 
between chains, respectively. The couplings considered 
are isotropic and antiferromagnetic, that is, Jii = J±_ = J 
and J > 0. As in the previous study ||, the ladders are 
investigated with a very efficient loop cluster algorithm 
Pj . The correlation length is obtained from the exponen- 
tial large-distance decay of the measured instantaneous 
spin-spin correlation function. We refine the original 
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simulations with open boundary conditions and re- 
peat them with periodic boundary conditions. These cal- 
culations yield low-T correlation lengths, £/a, of 7.1(1) 
and 10.3(1) for ladders of width n = 4 for periodic and 
open boundary conditions, respectively. The results for 
ri = 6 are 30.5(10) and 32.0(10) for periodic and open 
boundary conditions, respectively. Our corresponding 
analytical results for the periodic boundary condition are 
6.2 and 26.2 for n — 4 and n — 6, respectively |5|,0]. The 
spin gaps are obtained from independent numerical mea- 
surements of the uniform susceptibility. For n — 4 they 
are 0.234(4) and 0.160(4) for periodic and open bound- 
ary conditions, respectively. Similarly, for n — 6, we 
obtain 0.055(4) and 0.053(4), respectively, for periodic 
and open boundary conditions. It is seen that Eq. (Q) 
holds remarkably well. The temperature dependence of 
the correlation lengths obtained from Monte Carlo sim- 
ulations are shown in Fig. (0). For n = 6, the choice of 
boundary conditions indeed hardly matters. 
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FIG. 1. Correlation lengths obtained from Monte Carlo 
simulations as functions of temperature for spin-1/2 ladders 
composed of four and six chains. 

By studying the effective Euclidean action of the 
(2+l)-dimensional non-linear cr-model, 



f = ft L dT 1 lL ' 1 ''" 



(3) 



where /i is to be summed over the spatial directions, x 
and y, we calculate the temperature dependence of the 
correlation length in three distinct regimes characterized 
by the ratio /3hc/L y . The parameter is the bare spin 
stiffness constant at the spatial cutoff A -1 of the model, 
c is the spin-wave velocity defined on the same scale 
and fi, the staggered order parameter field, is a three- 
component unit vector. In the low-T (fihc ^> L y ) and 
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FIG. 2. The three different temperature regimes: a) low-T, 
b) intermediate-T, and c) high-T. The figures show the exten- 
sions of the two finite directions. For a) and c), the effective 
model is a (1+1) -dimensional quantum non-linear cr-model 
with fields constant in the shortest directions, and a coupling 
modified by the fluctuations in this direction. For b), the 
fields are constant along both the y— and r-directions, and 
the effective model is a one-dimensional classical cr-model. 

high-T (f3hc <C L y ) regimes we map (||) to a (1+1) 
dimensional quantum nonlinear cr-model where the spa- 
tial dimension refers to the a;-direction, and the extent of 
the "time-direction" is (3hc and L yi respectively. For the 
intermediate-T (fine ~ L y ) regime the effective model is 
the one-dimensional classical non-linear cr-model. The 
fluctuations that are integrated out are approximated by 
a one-loop modification of the coupling constant for the 
effective theory. Figure (||) illustrates the situations in 
the three different regimes. 

At T = the "Lorentz" -invariance of the non-linear 
cr-model gives the simple relation Eq. (|]) between the 
correlation length and the gap. For T ^ it is not so 
clear that the spin-wave velocity will survive quantum 
fluctuations unrenormalized. However, it is certainly the 
case that the quantum fluctuations at the one-loop level 
cannot renormalize the spin-wave velocity at any temper- 
ature. Therefore, within our one-loop approximation, it 
is consistent to promote Eq. (Q) to be valid at all tem- 
peratures, which we will do. 

For low temperatures, the same procedure as that used 
for the T = calculation ||, of integrating out the fluc- 
tuations along the finite-width (y-direction) to one-loop 
order, should work. However, at non-zero temperatures 
the resulting effective model is not the two-dimensional 
classical non-linear cr-model, but a (l+l)-dimensional 
quantum non-linear cr-model, 
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e(L v ,T) e(L y ,0) 
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exp phcA^/k 2 + *4 



(5) 
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The effects of the fluctuations in the y-direction are 
taken into account by modifying the coupling constant 
to one-loop order; for convenience, we separate out the 
temperature-independent part, e(L,0). In the above ex- 
pression k m = 27rm/LA, and the prime on the summa- 
tion sign means that the m — term is omitted from the 
sum. This expression is valid for ladders with periodic 
boundary conditions in all directions. In this Letter we 
only consider spin ladders with these boundary condi- 
tions. A short wavelength cutoff is imposed on the inte- 
gration and on the summation such that \Jk 2 + k 2 n < 1, 
where A is the cutoff scale chosen to conserve the area of 
the first Brillouin zone: Aa = 2^/tt. The results are only 
weakly dependent on this choice. Because of the absence 
of the topological term, the above action only describes 
spin-1/2 ladders of even width. 

The correlation length for the above model (^) can 
be obtained from a simple self-consistent approach. By 
relaxing the unitarity condition on f2 and adding a mass 
term, the condition < £1 2 >= 1 yields 



f i „ coth Jk 2 + A 2 (L y ,T)f3hcA/2 
Ne(L v ,T) f- V , '- = 1. 



k 2 + A 2 (L y ,T) 



(6) 



where N = 3 and A 
T = gives 

A(L,0) = 



(A£) 1 . Solving this equation at 



(7) 



sinh (2iT/Ne(L y ,0)) 



Unfortunately this is not the correct zero-temperature 
result (|J^]. To remedy this we can let Ne(L y ,0) be a 
free parameter and adjust it such that A(L y ,0) agrees 
with the correct zero-T result ||. Fixing Ne(L y ,0) in 
this way we can solve Eq. (6) numerically for the cor- 
relation length at non-zero temperatures. Due to the 
slight discrepancy between the T = correlation lengths 
given in Q and the correlation lengths obtained in the 
Monte Carlo simulations we adjust Ne(L y , 0) such that 
the T = correlation length agrees with the Monte Carlo 
data, that is £/a=7.1 and 30.5 for n = 4 and n = 6, 
respectively. The finite- T correlation lengths for these 
ladders are shown in region I in Fig. (^|) . 

For intermediate temperatures, (3%c r 
tuations in both the y- and the r-directions are inte- 
grated out. The low-energy effective action is then a 
one-dimensional classical non-linear cr-model, 
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where the coupling constant to one-loop order is 



t-\L y ,T) = PLpl- J g 
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FIG. 3. Correlation lengths as functions of temperature in 
the three distinct regimes for spin ladders of width n — 4 and 
n — 6, compared with the data from numerical simulations. 
The three regimes are labeled: I) low-T, II) intermediate-T, 
and III) high-T. 



Here, q m = 2irm/L y , uj n = 2im/ (3hc, and the prime on 
the sum indicates that the m = n = term is omitted. 
Using a combination of Ewald and Poisson summation 
techniques, as in Ref. @, we get 



t-\L y ,T) = ^ + L y A(^) 

4.7 L y 



(10) 



where £j is the Josephson length in the Nccl phase [[|, 
and the scaling function A is 



Z7T 



where 



dyy- 1 ' 2 {X{-Kxy)X{-Ky/x) - 1) - 2 



x (y)= E 



(ii) 



(12) 



The onc-dimensional classical cr-model is nothing but 
the continuum limit of the classical Heisenberg spin chain 
on a lattice, which has the action 



s = -Kj2& i -n i+1 . 



(13) 



The correlation length for this lattice model has been 
calculated exactly in pl[, and is 



1 



In [cothX -l/K] 



(14) 
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Taking the continuum limit a — > 0, with Ka kept fixed, 
we find 



5=—^— / dx[d x Cl) + const. 



(15) 



The one-dimensional non-linear cr-model is a finite the- 
ory for which no ultraviolet regularization is necessary, 
and there are no ambiguities associated with taking the 
continuum limit. Thus from Eqs. (^|) and ( [l5|) we can 
identify Ka = \/t(L y , T). Taking the continuum limit of 
Eq. (PJ), with Ka held fixed, we obtain 



e _ t-\L y ,T) 



(16) 



We note that this expression has no adjustable parame- 
ters. The results for ladders composed of four and six 
chains are shown in region II of Fig. (||). 

For temperatures such that L y 3> (3Tic, it is reasonable 
to map the action Eq. (|^) to a (l+l)-dimensional quan- 
tum non linear a -model in which the "time direction" is 
really the width L y — na of the ladders, 
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n 2 f (L v ,T) J dx ] Q dy 



(17) 



Here, the effective coupling constant is 
1 1 



f(L v ,T) /(oo,T) 



dk ^-v 
2^ ^ 



exp LyK^Jk 2 + uj 2 - 1 



(18) 



and w„ = 2Trn/ (3hcA. This is very similar to the low-T 
case; one just has to interchange f3hc with L y . In the 
low-temperature regime, we imposed a cutoff in the two 
spatial directions. Here, it is convenient to impose the 
cutoff in the x — r-plane. We have checked that this 
change in the cutoff procedure does not alter the results. 
In the same way as for the low-temperature regime we 
find the correlation length by a self-consistent equation, 



Nf(L y , T) 



i 1 

1 dk 



coth Jk 2 + A 2 (L y ,T)L y A/2 



2tt 



A 2 (L y ,T) 



(19) 



The quantity Nf(oo, T) is adjusted such that the solu- 
tion of the equation at L y = oo, £(oo,T), corresponds to 
the correlation length of the (2+l)-dimensional quantum 
antiferromagnet which was calculated in |? 13 1, 



f(oo,T) _ eV2Z c 



%TtSZ n 



^S 2 Z Ps J/T 



AttS 2 Z Ps J 



where Z c = 1.19 and Z Ps = 0.740 ]lj. Having fixed 
Nf(oo,T), we calculate k{L y ,T) by solving (|l|) numer- 
ically. The results for n = 4 and n = 6 correspond to 
region III of Fig. (§). 

It is evident from Fig. (^) that the analytic results 
for the correlation lengths agree remarkably well with 
the numerical simulations on even-legged spin-1/2 lad- 
ders with periodic boundary conditions. Since our effec- 
tive non-linear cr-model does not contain the topological 
term, our results are strictly valid for even-leg ladders. 
However, the topological term is important only for the 
physics H at very low energies. Thus, we expect that our 
expressions for the intermediate- and high-temperature 
regimes are valid also for odd-leg ladders. It would, how- 
ever, be interesting to have a more quantitative under- 
standing of the crossover to the WZW-model for odd- 
legged spin ladders. 
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